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Abstract. 

Entanglement of the excitonic states in the system of two coupled semiconductor 
microcrystallites, whose sizes are much larger than the Bohr radius of exciton in 
bulk semiconductor but smaller than the relevant optical wavelength, is quantified 
in terms of the entropy of entanglement. It is observed that the nonlinear 
interaction between excitons increases the maximum values of the entropy of the 
entanglement more than that of the linear coupling model. Therefore, a system of 
two coupled microcrystallites can be used as a good source of entanglement with 
fixed exciton number. The relationship between the entropy of the entanglement 
and the population imbalance of two microcrystallites is numerically shown and 
the uppermost envelope function for them is estimated by applying the Jaynes 
principle. 
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1. Introduction 

Recently, quantum computation and quantum information excite large enthusiasm 
among theoretical and experimental physicists. Many theoretical and experimental 
researches have been devoted to the preparation and measurement of the entangled 
states that have been considered to be a key ingredient in the realization of the 
quantum computation and information. 

We can quantitatively understand the role of quantum entanglement as a resource 
for communication by studying the quantum teleportation, which was first realized 
for discrete variables and then for continuous variables The possibility for 
the generation of the maximally entangled states with a fixed photon number from 
squeezed vacuum states or from mixed Gaussian continuous entangled states by the 
quantum non-demolition measurement [Sj has been theoretically studied. Quantum 
teleportation using an entangled source of fixed total photon number has also been 
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theoretically investigated in [3]. However, it has been shown that a generation of the 
maximally entangled states with a fixed finite photon number is a great challenge to 
modern technology. So the exploration of new maximally entangled bipartite source 
with fixed particle number is very interesting and important from both experimental 
and theoretical view points. 

Several schemes based on coupled quantum dots have been proposed for 
fabricating quantum gates [HE]. The entanglement of the exciton states in a single 
quantum dot or in a quantum dot molecule has been demonstrated experimentally 
[71 |S] . The authors of references I1UI 111! IT^) theoretically investigated the 
entanglement of excitonic states in the system of the optically driven coupled quantum 
dots and propose models to prepare Bell, Greenberger-Horne-Zeilinger (GHZ) and W 
states. The above investigations mainly focus on the case in which the sizes of quantum 
dots are smaller than the Bohr radius of exciton in bulk semiconductor, so that the 
quantum dots have discrete energy levels and one energy state can only be occupied 
by one exciton because of the Pauli principle. In such quantum dot systems, it is easy 
to define a qubit by different spins of excitons, or by single-exciton and no-exciton 
states. So when we want to encode quantum information in systems of quantum dots 
by a finite-dimensional quantum state space (called qudit system) , we need to consider 
higher energy levels. However, if we can enlarge the size R of the quantum dot so that it 
is larger than the Bohr radius as of the exciton in the bulk semiconductor but smaller 
than the relevant light wavelength A, in such a system of quantum microcrystallites 
(which can also be called large quantum dots system ^HIEHEI), the excitons can 
occupy the same level when the average number of excitons per Bohr radius volume 
is less than one. Then we can encode quantum information in the finite Hilbert 
space built by the different exciton number states. Quantum information properties 
of semiconductor microcrystallites coupled by a cavity field were investigated by us in 
Refs. |18| where, in particular, we have proposed a realization of symmetric sharing 
of entanglement for the exciton states in semiconductor microcrystallites and also 
studied the environment effect on the qubit. 

Motivated by these considerations, in this paper we will study a system of 
two semiconductor microcrystallites coupled by Coulomb interaction. We will study 
the entanglement of the qudit excitonic states in these two coupled semiconductor 
microcrystallites with fixed exciton number. 

The paper is organized as follows: In Sec. II, we first give a theoretical description 
of two coupled semiconductor microcrystallites. We obtain the analytical solution of 
the system by virtue of the Schwinger representation of the angular momentum. In 
Sec. Ill, the entanglement of the two subsystems is discussed using the von Neumann 
entropy under certain initial conditions. We discuss and numerically analyze the 
oscillation of the exciton population imbalance between two quantum dots. The 
relationship between the entanglement and oscillation of the imbalance population 
of excitons is demonstrated in detail in Sec. IV. Finally, we give our conclusions and 
some comments. 

2. The model and its analytical solution 

We consider a system which consists of two semiconductor microcrystallites coupled 
by Coulomb interaction. We assume that two semiconductor microcrystallites are 
completely symmetric. Their sizes are larger than the Bohr radius as of an exciton 
in a bulk semiconductor but smaller than the relevant light wavelength A. Using two- 
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band approximation to model these two coupled semiconductor microcrystallites, we 
have the Hamiltonian 

2 2 
H = H X]ka]a k + ^2 Xjki m a]a\aia m , (1) 

j,k—l j,k,l : m=l 

where a] (a 3 -) are creation (annihilation) operators of excitons, which are electron-hole 
pairs bound by the Coulomb interaction. We assume that the density of the excitons 
is so low and the external confinement potential to microcrystallite is so weak that 
the exciton operators Oj (aj) can be described approximately as bosonic operators, 

that is they satisfy the commutation relations of the ideal bosons, [ojjfflt] = Sjk> 
which is somewhat different from that in Ref. ^Hj- The label j=l (2) denotes the 
microcrystallite A (B). In the Hamiltonian QJ, the deviation of the exciton operators 
from the ideal bosonic model is corrected by introducing an effective nonlinear 
interaction between the hypothetical ideal bosons. In general, the parameters Xjk 
are different from each other, however in this study we consider two completely 
equivalent microcrystallites, which have nearly the same Bohr radius of the excitons 
and transitional dipole moment. For the sake of simplicity, we assume the parameters 
Xjk — lu > for j — k, which means that the two semiconductor microcrystallites 
have the same transition frequency from the valence band to the conduction band, 
and assume positive real numbers Xjk = 9 f° r 3 7^ which correspond to a linear 
coupling constant of the two semiconductor microcrystallites. The parameters Xjklm 
of the nonlinear terms are also taken as positive constant and set to be Xjklm = X- 
Under these assumptions, the Hamiltonian 0} can be simplified as follows 

H = hVLN + h X N 2 + hG(a\a 2 + 4«i) 

+ hx(a\a 2 + 4m) 2 , (2) 

where f] = uj — 2x, and G = g — 2x + 2%N. It is obvious that N = a\ai + a 2 a2 
is a constant of motion, which means [N, H] = and the total exciton number of 
the coupled semiconductor microcrystallites is conserved. In such a situation, we find 
that it is more convenient to get the solution of the Schrodinger equation governed 
by Hamiltonian J5J by virtue of Schwinger representation of the angular momentum 
J201 121 j . That is, we can introduce the angular momentum operators 

Jx = ^(a\a2 + 4ai), (3a) 
\-.{a\a 2 - a\a x ), (3b) 



Jy 



2i 



J* = ^( a l fl l _ a 2 fl 2) (3c) 

from the bosonic annihilation and creation operators of the two exciton modes. The 
operators of (|3ajl - l|3*cf satisfy the commutation relations for the Lie algebra of SU(2): 

[Jj,Jk]= iSjkiJu j,k,l = x,y,z, (4) 

where the Levi-Civita tensor Sju is equal to +1 and —1 for the even and the odd 
permutation of its indices, respectively, and zero otherwise. From l|3ajl-l|3c|l. the total 
angular momentum operator can be expressed as 

^ = y(y + l). (5) 
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In fact, N itself commutes with all the operators l|3af> - l|5c) l and (JSJ). For a fixed total 
excitonic number N ', the common eigenstates of J 2 and J z are the two-mode Fock 
states H3E] 

\j,m)s = K,m 2 ) = °! J P J O) (6) 

vu + m )'u - m )' 

with eigenvalues j = AT/2 and m = —Af/2,---,Af/2, where |mi,m2) is the Fock 
state with mi = j + m excitons and mi = j — m excitons in microcrystallite A and 
microcrystallite B respectively. Although j ± m (7711,7712) must be integers, j and m 
can both be integers or half-odd integers. For consistency, j is replaced by AT/2 in the 
following expressions. For clarity, the subscript s is used to indicate the Schwinger 
angular momentum basis in © and the following equations. By using equation we 
take j = 1 and m = as a detailed example to show the relation between the Schwinger 
angular momentum basis and Fock state basis. That is, |1, 0) s = aja^O) = |1, 1) which 
means that when the quantum number of the total angular momentum j is 1 and its 
z component m is 0, there is one exciton in each microcrystallite respectively. 

In terms of an SO(3) rotation e t7T ^ 2Jy of 2hQJ x + Ah\ J 2 , equation (0) can be 
simplified into : 

H = Mm + a X N 2 + 2%GJ X + Ah X J 2 

= h£lN + HxN 2 + 2hGe- l{ " T/2)J yJ z e ii7r/2)J y 

+ 4h X e- li7r/2)J yJ 2 e li7r/2)J y. (7) 

The eigenfunctions ^ m /2,m an d the eigenvalues Ej^/ 2 ,m °f Hamiltonian Q can be 
obtained easily as 

K 

\^K m ) s ^e^ 2 ^\^,m) s = J2 ^(DlyV)- (8a) 

Ess/2,m = m Af + h X Af 2 + 2hQm + 4h X m 2 , (8b) 
where AT denotes the total excitonic number of two microcrystallites; Q — g—2x+2\M , 
and Wigner's formula for f^/ m (f ) reads as 

^ m (?) = E(- 1 )^ m_m '(5) f 



yj(% + m)!(^ - m)!(^ + m')!(^ 



~ + m - fc)!fe!(=y - Jfe - m')!(fc - m + m')\ ' ^ 

where we take the sum over k so that none of the arguments of factorials in the 
denominator is negative. The wave function |\&(f)) of the system for the initial 
condition \^>(t = 0)) is given by 

oo m=H/2 

m)) = eM-itH/n)\m) = E E 

A/"=0 m=-Af/2 

x exp (-^t) (*# J*(0))|*v m > a . (10) 



The coefficients ($« m |\I/(0)) are rotating matrix elements that can be determined by 
Wigner's formula. Equation lfTU|) is basic and will be used in our further discussions. 
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In the following two sections, we will discuss the entanglement of two exciton modes 
and demonstrate the oscillations of population imbalance for excitons between two 
quantum microcrystallites. 

3. Entanglement of the excitonic states 

Quantum entanglement plays the key role in the quantum information and quantum 
computation. In general, for any pure composite state \iP)ab of a bipartite 
system whose state space is Ha <8> Hb, the entanglement can be measured by von 
Neumann's entropy of any reduced density operator pa = T^b(\iP)ab Ab(V'I) or 
Pb = TrA.(\ip )ab ab^), where the reduced density operator of system A is obtained 
by tracing out system B and that of system B by tracing out system A. In the 
context of quantum computation and information, in particular in studies of quantum 
entanglement, the reduced von Neumann entropy of a bipartite pure state \iJj)ab is 
usually referred to as the entropy of entanglement or simply the entanglement (see, 

e.g., EHE3E1ES]): 

e (p) = - Tr (PA hipA) = -Tr(/9 B In p B ) = A„ In A„ (11) 

a 

where A„ are the eigenvalues of either pa or ps- They form the (square of the) 
coefficients of the Schmidt decomposition of the bipartite pure state. That is, the 
bipartite pure state \tp) ab can be expressed by a set of biorthogonal vectors using the 
Schmidt decomposition as 

\i>) = \1P) A B = V^~n\a n ) A \Pn)B (12) 

n 

where we have chosen the phases of our basis states so that no phases appear in 
the coefficients A„ in the sum of |Q , {\a n ),n = ■ • •} and {\(3 n ),n — ()■••} are 
orthonormal states of the two subsystems A and B respectively. 

The total number of excitons in the whole system is fixed by the initial condition, 
such as L. Based on such a condition, the maximally entangled state of the system is 

|A/) = —L= (13) 

where \L — 1,1) represents the state with L — I excitons in microcrystallite A and 
I excitons in microcrystallite B. According to definition the entropy of the 

entanglement for the maximally entangled state (|13|) is ln(l + L). 

Experimentally, the easier way is to excite some number of excitons in one 
microcrystallite, or both of microcrystallites. Without loss of generality, we first 
assume that microcrystallite A is initially excited with L excitons, and no excitons 
initially exist in microcrystallite B. So the modes of the two microcrystallites are 
disentangled at the initial time t = 0. That is, the state of the whole system is a 
tensor product of the states of two subsystems A and B, i.e. \fy(t = 0)) = |L)yt|0)s 
with the Schwinger realization of this initial state as \^(t = 0)) = | U, U) s - The system 
and each mode of the system are in pure states and have zero entropies. 

It is well known that any pure state remains pure with the unitary evolution, but 
it is not true for each subsystem. With the evolution, the initially pure state of each 
subsystem can be transformed into mixed states. If the von Neumann entropy E(p), 
given by Eq. 111(1 . is increased, then the two subsystems become more entangled. 
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Thus, the degree of entanglement between the two subsystems in the system at any 
time can be described by (JTTJ. 

For the initial state |4, U) s of the system, the total wave function of the system 
can be obtained from (jTU)l as 

l'=L/2 E 

|*(t)>= Y, ex P(^*)^|,rl2'2 )s| **^' )s ' (U) 

l'=-L/2 

where the normalized coefficients (^>lj,\^, ^) s are determined by (|5a|) - (fTU|l . We will 
use equation (fH|l to discuss the degree of the entanglement of the two subsystems for 
the following several concrete examples. First, we consider that there is initially one 
exciton in the microcrystallite A, i.e., L = 1. In this case, we can obtain the wave 
function of the whole system from 114|) and @) with the initial state |i, i) s as 

\^)(t)} = co S (gt)\^) s -isin(gt)\^~) s 

= cos{gt)\l) A \0) B -isin{gt)\0) A \l) B , (15) 

where subscripts A and B denote dots 1 and 2, respectively, and we have omitted the 
global phase factor of time dependence, e~^ n+2x ^. The entropy of entanglement can 
be calculated as 

E {1) {t) = -cos 2 (.gi) ln[cos 2 (.gi)] - sin 2 ( ff i) ln[sin 2 (#<)]. (16) 

It is seen that the entropy of the entanglement periodically evolves with zero values 
at times gt = k^ where k is an integer. The entropy of the entanglement reaches 
its maximum values In 2 at times gt = (2k + 1)? with integer k, then the maximally 
entangled states of the system is 

II l\_4- P -i2gt\l -V. , 

,(17) 



|^(l) (t)) = \hl)s + e- t29t \h-^)s = \l) A \0) B +e-^\0} A \l) B 



V2 y/2 
If there are initially two excitons in microcrystallite A, then the wave function of the 
whole system with the initial condition \^(t = 0)) = |1, l) s becomes 

= ai \l, -l) s + a 2 \l, 0) s + aa|l, 1), 

= ai|0)^|2) B + a 2 |l)^|l) B + a 3 \2) A \0) B (18a) 

and 

ai = I [ e i2st + e -»2( 3 +4 X )t _ 2 ] ^ ( 18b ) 

^^^[e- 12 ^ 4 ^*-^ 2 ^], (18c) 

a 3 = i [e' i2 f* + e -»2(s+4 X )t + 2 ] (18d) 

where the global phase factor of time dependence e~ l2 ( n+2x - lt has also been neglected. 
We can obtain the entropy of the entanglement as 

3 

E^(t) = -Y\an\ 2 \n\a n \ 2 , (19) 

where a n determined by ijl8b(l - i|18d|) . There is no solution for the time t in which 
the system exactly evolves into a maximally entangled state such that the entropy 
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of entanglement is In 3. But numerically we find that the difference between the 
maximum value of the entropy of the entanglement and the entropy of the maximally 
entangled state is of the order 10 -5 . At this point, we say we can nearly obtain the 
maximally entangled state when the total number of the excitons in the system is equal 
to two. Figures ^a-ndOl demonstrate the evolution of the entropy of the entanglement 
as a function of gt with two different sets of parameters for given different initial 
states. They show that when the second-order coupling between excitons is weak, the 
period for reaching the maximum values of the entropy of the entanglement becomes 
longer. It is worth pointing out that if the second-order coupling constant is equal to 
zero and only linear coupling term remains, then the exciton density is lower and the 
entanglement is diminished (see solid curve with dots in comparison to other curves 
in figure 3). So, in order to obtain high entanglement of excitons, we should increase 
the exciton density in the coupled microcrystallites. The above discussion also shows 
that the entanglement between two microcrystallites depends on both the ratio xl 9 
and initial conditions. By numerical calculations, we can find that the survival time 
of higher entanglement also depends on xl 9 and initial conditions. The survival 
time of higher entanglement is longer for the larger ratio xl 9 when there are two 
excitons initially excited in one of the microcrystallites, but it almost the same for 
different ratios xl 9 when there is one exciton initially excited in each microcrystallite, 
respectively. 

For an arbitrary number of excitons L, the wave function of the whole system is 
described by lJTI)l. Then the entropy of the entanglement can be calculated as 

L_ 

2 

E (L \t) = E^°\t) = - l/M 2 m|/M 2 (20a) 



with 



m ' = — 



Pm'= £ ^!, m (f)e^ £ *""^i, m (|), (20b) 



m= — 



where V^, m (§) is determined by 10. 

Equations (|20all20b|) are confined to the case of the excitons initially excited in 
one microcrystallite only. Experimentally, we can also initially excite excitons in both 
microcrystallites. If P and Q excitons are initially excited in microcrystallites A and 
B then we have j — (P + Q)/2 and m — (P — Q)/2, and so the initial state can be 
expressed as |<&(£ = 0)) = \P)a\Q)b = l^^j ^f^)a- At time t, the wave function of 
the whole system can be written as 

(P+Q)/2 

!*(*)>= E exp(-iW^) 
l'=-(P+Q)/2 

P 4- O P — o 

x (tfpja^l-p, — 2^).|* £ ja, I ,).. (21) 

We can also obtain the entropy of entanglement for the above initial condition as 

(P+Q)/z 

E ^Q){t) = - Wn>\ 2 M&>? (22a) 

n'=-(F+Q)/2 
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Figure 1. 

£(2,0) 

is plotted as a function of gt for (a) \l 9 = 0.34, (b) xl g = 0.01 
when two excitons are initially excited in one of the microcrystallites, and E^' 1 ' 
is plotted as a function of gt for (c) x/<? = 0.34, (d) \l 9 = 0.01 when there is one 
exciton initially in each microcrystallite. 



with 

(P+Q)/2 

" p+q _i£ p+(i p+Q n 

P'n> = E V 4^J^ e (f )• (22b) 

n=-(P+Q)/2 2 

Equations 122all22b|) for the entanglement were obtained from the wave function (|21ll 
being a complete solution of the model. Thus, the contributions of all elements (also 
those off-diagonal) of the density matrix are included. In figure the entropies 
of the entanglement E^ p ^ are depicted as a function of gt for different exciton 
numbers P and Q in the first and second microcrystallite, respectively, but for 
the same total number of five excitons in both microcrystallites. Numerical results 
show that the evolution behavior of the entropy of the entanglement for L = 5 is 
similar to that for L = 2. However, when we compare the maximum value of the 
entropy of entanglement with the entropy of the maximally entangled state, we see 
that a time t at which the generated entangled state is very close to a maximally 
entangled state cannot be found for any xl 9- The survival time of the maximum 
entanglement is different for different initial conditions and xl 9- The smaller difference 
of initial populations corresponds to longer time in a larger entanglement area for 
fixed xl 9 when the total number of excitons in two microcrystallites is five. We also 
find that the maximum entanglement for fixed xl 9 is the same for different initial 
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Figure 2. = E 1 - 5 '^, fit 4 ' 1 ) , and _B< 3 ' 2 > are plotted as a function of gt for 

x/g = 0.34 (left panel) and \l 9 = 0-01 (right panel). 



conditions for the total number L of excitons in two microcrystallites L = 2, 3 only, 

ie £(2,0) _ £,(1,1) _ £(0,2) _ , „ , ^(3,0) _ £(2,1) _ £(1,2) _ £(0,3) _ „, „ 
i.e. J^max — ^max — ^max — in o J-^max — ^max — ^max — ^max — aiixa. 

However, if L > 3 then £max™'™' ^ £mai m ' m ' for n ^= m < L/2. Nevertheless, the 
differences are relatively small - only at the second digit after the decimal point, as we 
have checked numerically up to ten excitons. Concluding, we find that the maximum 
values of the entropy of entanglement depends both on the initial conditions and on 
X/g when L > 3. 

We can use equations H20aH20b|) as an example to further discuss the entanglement 
of the excitons between the two subsystems A and B of the system for any initially 
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Figure 3. The maximal values E^^x of the time-dependent von Neumann 
entropy E^(t) are plotted as a function of exciton number L < 10 for xl 9 = 
(dots in solid curve), x/g = 0.01, 0.34, and 0.8 (squares in solid curve). The 
values of the von Neumann entropy of the maximally entangled states are marked 
by circles in the dotted curve for different exciton numbers L. 

given exciton number. Here the entropy of the entanglement of two subsystems is 
discussed up to ten excitons. An exact time i, which corresponds to a maximally 
entangled state, can be obtained when the system has one exciton. Numerically up to 
a evolution time of gt < 600, when the total number of excitons in our system is two 
or three, we find that the difference between the maximal entropy of entanglement 
and the entropy of the maximally entangled state is of the order 10 -5 for xl 9 ^ 0- So, 
the maximally entangled states can be approximately obtained. The maximum values 
^max of the entropy of the entanglement for the two subsystems are also numerically 
calculated by using H20al) - (|20b() up to the time gt < 600. Figure [3] shows a comparison 
of the maximal values E^stx of the entropy of the entanglement with the values of the 
entropy of the entanglement for maximally entangled states up to ten excitons for 
x/g = 0, 0.01, 0.34, 0.8. Figure [3| shows that the entropies of the entanglement for 
x/g = are smaller than those for xl 9 = 0.01, 0.34, 0.8 for all exciton numbers, 
which means our system can achieve a larger entanglement than that of the linear 
coupling model. At this point, this coupled microcrystallites system can be considered 
as a better source of entanglement with a fixed exciton number than that of the 
linear coupling model with the same particle number. For a fixed exciton number, we 
find that the difference in the maximum values of the entropies of the entanglement 
between any two different parameter ratios among xl 9 = 0.01, 0.34, 0.8 is very small. 
Every square in the solid curve of figure|31actually denotes the maximum values of the 
entropies of the entanglement for three different values of xl 9- By comparing plots for 
a small (x/g = 0.01) and a higher (0.34) values of the nonlinearity in figures^|2 and 
13 it is clearly seen that the global maxima of entanglement are practically independent 
of the ratio of xl 9 (except x/g), but dependent of the initial conditions. However, 
higher nonlinearity corresponding to a larger xl 9 enables generation of globally higher 
entanglement for shorter evolution times. 

4. Oscillations of exciton imbalance 

In the previous section, we show that the excitonic state of two microcrystallites 
system has different degree of entanglement with the evolution. The population of the 
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excitons in the microcrystallites A and B is not at equilibrium. The nano-technology 
opens up a new direction of research for direct measurements of exciton dynamics. 
Probing one exciton at a time has been demonstrated by using the nonlinear nano- 
optics 26 . It is possible for experimentalists to observe the population of the excitons 
in the microcrystallites. Considering that experimentally, it is easier to initially excite 
excitons in one microcrystallite, in this section our discussions are the case when 
there are L excitons excited initially in one of the microcrystallites. We will show that 
the excitonic population imbalance between two coupled microcrystallites exhibits 
a collapse and revival phenomenon. The time-dependent population difference of 
excitons between two microcrystallites can be written as 

AA^(i) = N[ L \t) - N { 2 L \t), (23) 

where the average exciton number of each microcr ystallite JV/ i) (*) = (a]ai). The 
evolution of the population difference of the excitons is similar for different initial 
conditions. We can discuss the population imbalance of the two microcrystallites 
using equation (|10|) , but here we only consider the initial condition at which there are 
L excitons in microcrystallite A and no exciton in microcrystallite B. From (|14f> . we 
can calculate the number of excitons in each microcrystallite as 

Af 3 (*) = f-(-iy £ m '\^\\ (24) 

ro'=— I 

where I — 1 represents the microcrystallite A and I = 2 denotes microcrystallite B; 
time-dependent functions (3 m i are given by equation l|20b(l . Then we can obtain the 
population difference of the excitons between the two microcrystallites as 

t 

AN^(t) = 2 £ ™'l/9™'| 2 - (25) 

When only one exciton is initially excited in the microcrystallite A, the exciton 
population difference is 

AN w (t) = cos(2 5 i), (26) 

which is a simple sinusoidal oscillation. The population difference of the exciton is zero 
at times gt = (2k + l)f , and AN^ 1 ' reaches a maximum at times gt = fc^, where k is 
an integer. Comparing this result to equation (|16l) . we know that the two subsystems 
reach the maximal entanglement when the exciton population difference of the two 
microcrystallites is zero, but when the entropy of the entanglement of the system is 
minimum, the population imbalance is maximum. 

When there are two excitons initially excited in microcrystallite A, we obtain 

AN {2) (t) = 2 cos(4 X t) cos[2gt + A X t] . (27) 

Equation 127|l indicates that the population difference periodically oscillates with 
time. We numerically give the evolution of the population difference as a function 
of gt in figure Figure 0] shows that the exciton population difference of the 
two microcrystallites exhibits the beat effect. If the second interaction between the 
excitons becomes weak, the period of the envelope becomes longer than that of the 
stronger second interaction between excitons. 
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Figure 4. The evolution of the population difference between two 

microcrystallites is shown as a function of gt for (a) \/g = 0.34 and (b) \j g = 0.01. 



(a) (b) 




1 2 3 4 5 20 40 60 



gt gt 

Figure 5. The evolution of the population difference AN^ between two 
microcrystallites is shown as a function of gt for (a) \/g = 0.34, and (b) 

x /g = 0.01. 



The general formula for the exciton population difference of two microcrystallites 
can be given as 

AN (L) (t) = Lcos L - 1 (4:xt)cos[2gt + 4:(L~ l) X t\. (28) 

With increasing exciton number, the numerical results show that exciton population 
difference evolution in time resembles a collapse and revival phenomenon. As an 
example, the evolution of the exciton population difference in a system of the coupled 
microcrystallites with a total number of L = 5 excitons is plotted in figure as a 
function of gt. The smaller is the ratio between the second coupling constant and the 
first coupling constant, the clearer is this phenomenon. 

Figure [H] shows a relationship between the entropy of the entanglement in our 
system and the exciton-population difference. To estimate the uppermost envelope 
functions of the curves in figure we apply the Jaynes principle of maximum 
entropy |27) . In the Jaynes formalism, the entropy can exclusively be expressed in 
terms of the mean values. The method corresponds to averaging over the generalized 
grand canonical ensemble of states satisfying given constraints. Thus, the Jaynes 
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formalism of calculating the maximum entropy strongly resembles the standard 
problems of entropy maximization in statistical mechanics and thermodynamics. It 
is worth noting that, although the maximum entropy principle "has nothing to do 
with the laws of physics" [27|, it is a useful and widely applied tool in physics, e.g., in 
the quantum state reconstruction from data obtained in a measurement of a physical 
process |28j . By applying the Jaynes principle we find the maximum entropy to be 
expressed as a function of the population difference as 



E l j L) (AN^) = -^2p n lnp n = -(L + AN^^X^+lnZ^, (29) 

where p n = exp(— nM L ')/Z^ L ' = x n /Z^ L ' are given in terms of Lagrange multiplier 
\( L ' or x — exp(— \( L ') and the generalized partition function 



Z (L) = X i 

X — 1 

The Lagrange multiplier A^ L ) can be calculated from 



AN^=f(x)=L + T ^--2-±t^. (31) 

To obtain an explicit formula for Ej L ^ one has to find the relation inverse to 
= f(x). By closer inspection of (|31|l . we conclude that the analytical solutions 
of x = f~ 1 (AN^) exist only for L=2, 3, and 4. For example, the desired solution of 
x = /^(AN^) f or N = 2 is given by 



x - exp(-AC*>) = AiV(2)+ / 16 -? 2 ( ) AiV(2)) -. (32) 
PV ; 2(2 -AM 2 )) v ; 

For brevity, we do not present our analytical, but rather complicated, solutions for 
L = 3 and 4. The solutions for L > 5 are to be found numerically. The curves marked 
by circles in figure 6 correspond to our analytical (L = 2, 3, 4) and numerical (L = 5) 
estimations of the maximum entropy of entanglement as a function of AN^ L \ It is 
seen, in agreement with our analysis of figure 3, that with the increasing number L of 
excitons, the time- maximized entanglement E^^(t), which can be generated in our 
system, increases in comparison to the lower-exciton systems Em&x(t) > E^^^) 
for the same population difference A7V = AiV (i) = AiV (L_1) ^ ± L (with any L) 
for a given n/g. The states described by Ej(AN^) for AiV = correspond to 
the maximally entangled states. The numerical results depicted in figure show that 
the maximum entropy of the entanglement always corresponds to the points of the 
exciton-population balance. However, the condition of the exciton-population balance 
does not always imply the maximum entropy of the entanglement. For any population- 
imbalance points, the entropy of the entanglement is smaller, and when the population 
difference reaches its maximum, the entropy of the entanglement vanishes. 
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AA/W AN< 5 > 

Figure 6. The entropy of entanglement is plotted as a function of the population 
difference of microcrystallites for x/d = 0-34 when the total number L of excitons 
is 2, 3, 4, or 5. The uppermost circles represent values estimated by the Jayncs 
principle. 



5. Conclusions 

We have discussed the entanglement of the excitonic states in the system of the coupled 
microcrystallites with a fixed total exciton number by the entropy of the entanglement. 
When the total number of excitons is one, the maximally entangled state can exactly 
be obtained. The maximally entangled states can approximately be generated for two 
or three excitons. However, when the exciton number is more than three, we cannot 
obtain the maximally entangled state, but it is observed that the nonlinear interaction 
between excitons increases the maximum values of the entropy of the entanglement 
to higher values than those of the other models based on the linear coupling. So this 
exciton system can be used as a good source of entanglement. We also find that the 
entanglement between the two microcrystallites depends on the initial conditions, the 
different initial conditions will result in different entanglements. 

The oscillations of exciton population between two microcrystallites have also 
been investigated analytically and numerically If the system is prepared with 
one exciton, the population imbalance between two microcrystallites reaches the 
maximum when the entanglement of the two microcrystallites is the minimum; and 
the population reaches the balance when the entropy of the entanglement of the 
system is maximum. But when the exciton number is more than two, the numerical 
results showed that we could find the maximum entropy of the entanglement at the 
points of the exciton population balance. At some larger population imbalance points, 
the entropy of the entanglement was found to be smaller. The reason is that when 
the population imbalance becomes the largest, all excitons are concentrated on one 
microcrystallite, then the entanglement between two microcrystallites decreases. But 
in the collapse area, the average population of each microcrystallite is equal, so it is 
possible to find the maximal entropy of the entanglement. We hope these phenomena 
can be observed by experimentalists with the development of the nano-technology. 
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The main conclusion of our paper is that nonlinear interactions, resulting from 
high exciton density between microcrystallites, can increase exciton entanglement in 
comparison to a model of linearly-coupled microcrystallites with lower exciton density. 
We believe that this result could be important for semiconductor-based quantum 
information processing. Our conclusion is not only valid for the studied system, but 
can also be generalized to other Kerr nonlinear vs linear interaction models. 
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